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BACKGROUND 


Many  modem  explosive  applications  require  finite  element  and  finite  difference 
modeling  of  both  overdriven  detonation  and  lower  pressure  detonation  products  expansion. 
Large  detonation  wave  shaping  and  multiple  point  initiation  are  typical  overdriven  detonation 
applications.  Lower  pressure  detonation  products  expansion  is  required  for  material 
acceleration  applications,  such  as  explosively  formed  penetrators  and  shaped  charges.  Current 
thermodynamic  equation  of  states  used  in  dynamic  finite  element  and  finite  difference 
programs  are  either  parameterized  to  give  agreement  with  thermochemical  calculations  [1]  or 
experimental  copper  cylinder  explosive  expansion  experiments  [2,3].  Thermochemical 
calculations  have  proven  to  be  very  useful  for  the  prediction  of  explosive  products  properties, 
particularly  near  and  above  the  Chapman-Jouguet  state  [4,5].  Unfortunately,  they  do  not 
reproduce  the  products  expansion  behavior  accurately  enough  for  typical  warheads  design. 
Currently,  thermodynamic  equations  of  state  (JWL,  Wilkens)  [2,3]  used  for  warheads  design 
are  normally  calibrated  to  give  agreement  with  copper  cylinder  explosive  expansion 
experiments.  These  equations  of  state  have  not  been  calibrated  for  high  pressures  above  the 
Chapman-Jouguet  state.  Experimentation  [6]  and  comparison  with  thermochemical 
calculations  (Figures  1,2,  and  3)  have  demonstrated  that  a  poor  description  of  the  high  pressure 
region  exists.  In  order  to  achieve  a  suitable  equation  of  state  for  both  overdriven  and  lower 
pressure  products  expansion,  an  appropriate  equation  of  state  form  has  been  derived.  A 
previous  preliminary  report  [7]  briefly  describes  the  equation  of  state  and  parameterization 
methodology,  but  does  not  include  details.  This  report  provides  a  detailed  theoretical 
background  and  equation  of  state  formulation.  A  forthcoming  report  will  provide  a  detailed 
description  of  the  parameterization  methodology. 

EQUATION  OF  STATE  FORMULATION 


The  equation  of  state  form  was  chosen  so  as  to  adequately  describe  the  high  pressure 
regime  produced  by  overdriven  detonation,  and  yet  retain  the  low  pressure  expansion  behavior 
required  for  standard  material  acceleration  modeling.  To  this  end,  the  derived  form  is  based 
on  the  Jones-Wilkens-Lee  (JWL)  equation  of  state  [1]  due  to  its  computational  robustness,  and 
asymptotic  approach  to  an  ideal  gas  at  high  expansions.  Additional  exponential  terms  and  a 
variable  Gruneisen  parameter  have  been  added  to  adequately  describe  the  high  pressure  region 
above  the  Chapman-Jouguet  state.  The  resulting  equation  of  state  form,  named 
Jones-Wilkens-Lee-Baker  (JWLB),  is 
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For  consistency  with  the  JWL  equation  of  state,  V  is  defined  as  a  specific  volume  ratio,  V  a 


pb/p  and  E  is  defined  as  Em  poe  where  e  is  the  specific  internal  energy.  The  JWLB  equation 
of  state  form  is  based  on  a  first  order  expansion  around  the  principal  isentrope: 
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Using  the  Gruneisen  Parameter, 


Some  important  characteristics  of  the  equation  of  state  are  that  the  Gruneisen  Parameter  A,  is 
represented  as  an  analytic  function  of  specific  volume,  V.  A  +  1  approaches  a  constant 
V  dP 

adiabatic  gamma,  p  =  a»  1  for  large  V,  so  that  ideal  gas  behavior  is  asymptotically 

s 

approached.  The  principal  isentrope  description  is  essentially  identical  to  JWL,  with  the 
exception  of  an  increased  number  of  exponential  terms.  It  has  been  found  that  for  most 
explosives,  three  exponential  terms  (instead  of  two  in  JWL)  are  adequate  to  describe  the 
principal  isentrope  over  both  the  high  pressure  region  above  the  Chapman-Jouguet  state  and 
the  lower  pressure  expansion  region.  It  is  important  to  note  that  the  internal  energy 
referencing  is  defined  by  (7).  The  value  of  Eo  must  be.  consistent  so  that, 
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Eo  =  ^  e  RiVcj  +  §  VcT  -  i(Pcj  +  P0)(V0  -  Vcj)  . 

i  1 

Normally,  the  equation  of  state  parameters  are  chosen  so  that  Eo  has  the  value  Eo  =  pt,A H, 
where  A H  is  the  heat  of  detonation.  This  is  consistent  with  the  initial  internal  energy  of  the 
unreacted  material  having  a  value  of  zero. 

SPEED  OF  SOUND 

The  generalized  speed  of  sound  is  often  required  for  implementation  into  a  finite  element 
or  finite  difference  program.  The  speed  of  sound  may  be  easily  derived  as  follows, 
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ADIABATIC  GAMMA 


Another  useful  quantity  is  the  generalized  adiabatic  gamma.  The  adiabatic  gamma  may  also 
be  easily  derived, 
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PRINaPAL  ISENTROPE  PROPERTIES 

Often,  properties  along  the  isentrope  that  passes  through  the  Chapman-Jouguet  state  are  of 
particular  interest.  The  Gruneisen  parameter  is  given  by  (2).  The  adiabatic  gamma  is  given 
vj< 
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(3),  (15)  =>  ys  = 
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REACPED  PRODUCrS  HUGONiar 

Another  often  used  state  space  locus  is  the  reacted  products  Hugoniot.  Assuming  the  initial 
pressure  to  be  zero,  conservation  gives, 


Mass:  fxD  =  p(D-u) 
Momentum:  P  =  pJXu 
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CONCLUSION 

An  advanccvi  thermodynamic  equation  of  state  applicable  to  problems  involving 
overdriven  detonation  and  material  acceleration  has  been  developed  and  calibrated  for  several 
explosives.  This  manuscript  has  presented  a  theoretical  background  and  formulation  for  the 
new  equation  of  state.  A  forthcoming  report  will  describe  the  parameterization  methodology. 
The  new  equation  of  stme  maintains  required  low  pressure  expansion  behavior  while  providing 
a  better  high  pressure  description  applicable  to  overdriven  detonation.  Typical  applications  of 
the  new  equation  of  state  include  wave  shaped  and  peripherally  initiated  munitions.  Finally, 
the  equation  of  state  asymptotically  approaches  the  constant  gamma  equation  of  state  at  high 
expansions.  It  is  therefore  believed  that  the  equation  of  state  will  be  applicable  to  extremely 
large  volume  expansion  applications,  such  as  blast,  without  a  loss  of  accuracy  in  the  higher 
pressure  regions.  In  order  to  use  the  equation  of  state  for  very  large  expansions,  proper 
calibration  to  very  large  volume  expansion  experimentation  will  be  required. 


SYMBOLS 


P  a  Pressure 
p  s  Density 

V  a  Ratio  of  Specific  Volume  to  Initial  Specific 
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Volume 


E  m  Specific  Interna!  Energy  divided  by  Initial 
Specific  Volume 

X  =  Gruneisen  Parameter 

y  f  Adiabatic  Gamma 

c  =  Sound  Speed 

D  =  Detonation  Velocity 

u  s  Mass  Velocity 

A U  s  Pleat  of  Detonation 


Ap  Rp  C>  AAi’ 
BAi’  RAi’  10  ^ 


Equation  of  State  Constants 


%  s  Isentropic 

cj  s  Chapman-Jouguet  State 

u  ~  Initial  Conditions 
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EOS  Calculations 
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Figure  1.  Pressure  versus  specific  volume  for  the  principal  isentrope  of  octol  75/25  below  the 
Chapman-Jouguet  state.  The  thermochemical  calculations  (13KWK  and  JCZd)  agree  fairly  well 
with  the  standard  JWL. 


EOS  Calculations 


Figure  2.  Pressure  versus  specific  volume  for  the  principal  isentrope  and  reactive  Hugoniot  ot 
octol  75/25  above  the  Chapman-Jouguet  state.  The  standard  JWL  underpredicts  the  high 
pressure  region. 
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Figure  3.  The  Gruneisen  parameter  versus  specific  volume  for  the  principal  isemropc  and 
reactive  Hugoniot. 
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